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Convective and radiative heating environments have been computed for a three-dimen- 
sional ellipsled configuration which would perform an aerocapture maneuver at Neptune. 

This work was performed as part of a one-year Neptune aerocapture spacecraft systems 
study that also included analyses of trajectories, atmospheric modeling, aerodynamics, struc- 
tural design, and other disciplines. Complementary heating analyses were conducted by sep- 
arate teams using independent sets of aero thermodynamic modeling tools (i.e. Navier-Stokes 
and radiation transport codes). Environments were generated for a large 5.50 m length 
ellipsled and a small 2.88 m length ellipsled. Radiative heating was found to contribute up to 
80% of the total heating rate at the ellipsled nose depending on the trajectory point. Good 
agreement between convective heating predictions from the two Navier-Stokes solvers was 
obtained. However, the radiation analysis revealed several uncertainties in the computa- 
tional models employed in both sets of codes, as well as large differences between the pre- 
dicted radiative heating rates. 

p = density (kg/m 3 ) 

Subscripts 

rad = radiative 

conv = convective 

INTRODUCTION AND BACKGROUND 

A one year, multi-disciplinary study of a mission 
to Neptune in which aerocapture would be used to 
decelerate into orbit has been conducted. Computa- 
tional analyses of the convective and radiative aero- 
heating environments which the vehicle would 
experience are detailed herein, and results from other 
disciplines are presented in several companion 
papers 1 " 5 . 

Neptune Aerocapture Mission Concept 

In an aerocapture mission, atmospheric drag is 
employed in place of a conventional propulsion system 
to decelerate the vehicle into orbit (Fig. 1). Aerocap- 


NOMENCLATURE 

= reference area (m 2 ) 

= coefficient in reaction rate equation 
= drag coefficient 

= forward reaction rate (cm 3 /mole/s) 

= lift-to-drag ratio 
= mass (kg) 

= exponent in reaction rate equation 
= heat-transfer rate (W/cm ) 

= time (s) 

= temperature (K) 

= reaction temperature (K) 

= vibrational temperature (K) 

= axial distance measured from nose (m) 
= angle-of-attack (deg) 

= ballistic coefficient (kg/m ) 

(3 = m/(C D A) 

= activation temperature (K) 
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ture can result in large mass savings in comparison to 
propulsive deceleration. For this study a reference mis- 
sion concept 1 was developed for a 2017 launch with a 
10 year transit to Neptune of an orbiter designed for a 
scientific investigation of Neptune and its moon Triton. 
For the reference mission guidelines, it was determined 
that the mass savings resulting from an aerocapture 
maneuver would be necessary to deliver the required 
payload. 


Atmospheric Entry Interface 



Figure 1: Illustration of Aerocapture Mission 
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Vehicle Configuration 

An aerodynamic trade-off study 2 was conducted 
to define the shape of the vehicle’s aeroshell. The 
design objectives were to achieve a lift-to-drag (L/D) 
ratio of 0.8, minimize the ballistic coefficient ((3), max- 
imize the volumetric efficiency, and fit within the 
launch vehicle shroud. The configuration selected was 
the “flattened ellipsled” geometry shown in Fig. 2. A 
basic ellipsled configuration can be defined by an ellip- 
soid nose section followed by an elliptical cross-section 
cylinder. The basic ellipsled can then be “flattened” by 
shrinking the minor axis of the bottom half of the ellip- 
tical cross-section of the vehicle. 

The dimensions shown in Fig. 2 define the geome- 
try of a 5.50 m long vehicle which was the focus of the 
first phase of this study. Results from this phase of the 
study were used to conduct a design iteration, which 
resulted in a new, smaller vehicle. This second geome- 
try was a scaling of 52.36% from the original design 
which produced a 2.88 m long vehicle. 

Vehicle Trajectories 

Atmospheric trajectories for both vehicles were 
generated 3 using the Program for Optimization of Sim- 
ulated Trajectories 6 with a Neptune atmospheric 
model 4 . For the large vehicle, a ballistic coefficient ((3) 
of 400 kg/m 2 was used, while for the smaller vehicle a 
ballistic coefficient of 895 kg/m 2 was used. Reference 
convective heating rates were computed along these 
trajectories for a 1.00 m hemisphere using a Sut- 
ton-Graves 7 formulation. Several thousand trajectories 
were simulated, from which worst-cases for convective 
heat loads and heat rates were identified. These trajec- 
tories are shown in Figs. 3 and 4. For the original large 
vehicle, aeroheating predictions were generated at five 
points (including peak heating) along the max heat-rate 
trajectory, while for the final smaller vehicle heating 
rates were computed at the peak heating point along the 
max heat-rate trajectory and at seven points along the 
max heat-load trajectory. These heating predictions 
were then used to develop the aeroheating environ- 
ments required to design 5 a Thermal Protection System 
(TPS) for the vehicle. Free stream conditions for these 
points are given in Tables 1 and 2. The angle-of-attack 
for all cases was a = 40-deg. 
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Figure 3: Max Convective Heat-Rate Trajectory for 
Large Ellipsled 



t (sec) 

Figure 4: Max Convective Heat-Rate and Max 
Heat-Load Trajectories for Small Ellipsled 

COMPUTATIONAL TOOLS 

Convective-radiative heating environments were 
generated by two independent teams using separate 
flow field and radiation transport codes. Comparisons 
of results from the two analyses were performed for 
each set of trajectory computations in order to verify 
the results and identify sources of uncertainty. 

Within the context of this one-year study, the focus 
of the work was on the generation of environments for 
the system study rather than computational tool devel- 
opment. However, computation of the Neptune aero- 
capture environment posed several challenges which 
will be discussed in later sections, and further research 
into computational tools and models for the Neptune 
missions will be required to address these issues. 

Flow Field Solvers 

Flow field solutions were generated using the 
LAURA (Langley Aerothermodynamic Upwind Relax- 
ation Algorithm) 8,9 and DPLR 10 (Data Parallel Line 
Relaxation) codes. Both codes are three-dimensional, 


structured, finite-volume Navier-Stokes solvers and 
support multiple-block computations on distributed 
nodes using the Message Passing Interface (MPI). 
Inviscid fluxes are computed in LAURA using the Roe 
flux splitting 11 method with Harten’s entropy fix 12 and 
Yee’s symmetric total- variation diminishing limiting 13 , 
while a modified Steger- Warming flux vector 
splitting 14 with MUSCL extrapolation to third order 
with a minmod limiter 15 is implemented in DPLR. Pre- 
vious studies have shown that both codes produce simi- 
lar results when the same kinetic and transport 
properties models are implemented in each code (e.g. 
Refs. 16-17). 

A 13-species (H 2 , H, H+, He, He+ CH 4 , CH 3 , 
CH 2 , CH, C 2 , C, C+, e-) Neptune atmospheric model 
was employed for flow field computations with default 
free stream mass fractions of 0.6246 for H 2 , 0.2909 for 
He, and 0.0846 for CH 4 . A simpler 5-species model 
(H 2 , H, H+, He, e-) with freestream mass fractions of 
0.668 for H 2 and 0.332 for He was also employed for 
some numerical studies. Two reaction sets were used: 
the first was taken from Nelson 18 (with the addition of 
certain reactions from Park 19 and Leibowitz 20 ); in the 
second set, hydrogen and helium dissociation and ion- 
ization reactions were replaced with those from 
Leibowtiz 20 . Equations for the forward reaction rates 
for both sets are listed in Tables 3 and 4. Reverse rates 
were computed from the definition of the equilibrium 
constant, which was determined by evaluating the 
Gibbs free energy from thermodynamic data supplied 
by McBride 21 . 

A radiative wall equilibrium temperature boundary 
condition with a surface emissivity of 0.90 was 
imposed at the surface. “Super-catalytic” behavior 
(recombination to free stream mass fractions) was 
imposed to provide conservative heating estimates. 

Radiation Transport Codes 

Uncoupled radiative transport computations were 
performed using the flow field solutions as input data. 
An updated version of the RADEQUIL code 22,23 was 
used to process LAURA inputs and the NEQAIR96 24 
code was used to process DPLR inputs. 

RADEQUIL and NEQAIR96 are both used to cal- 
culate radiative emission and absorption from input 
flow field properties with the one-dimensional tan- 
gent- slab approximation in which it is assumed that 
radiation transport takes place only in the direction per- 
pendicular to the surface. The populations of the 
excited states of the various species are assumed to fol- 
low Boltzmann distributions. In RADEQUIL, a line 
group approximation is used to model atomic line tran- 
sitions and a smeared band approximation is used to 
model molecular transitions, whereas in NEQAIR96, 
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line-by-line computations are performed for all atomic 
and molecular transitions. In this study, transitions in 
hydrogen atoms due to excitation, bound-free 
photo-ionization, and free-free transition were consid- 
ered in both codes. Molecular transitions of H 2 were 
also included in the RADEQUIL analysis. The transi- 
tions occurring in C, C+, and C 2 were modeled, but 
were found to have negligible contributions due to the 
low concentrations of these species. 

RESULTS AND ANALYSIS 

Heating environments computed for the large and 
small ellipsled configurations are presented in this sec- 
tion. Uncertainties due to flow field and radiation 
transport solver implementation, vibrational non-equi- 
librium, atmospheric composition, radiation-flow field 
coupling, and radiation models are also discussed. 

5.50 m, P = 400 kg/m 2 Vehicle 

Convective heat transfer computations were per- 
formed for five points (Table 1) along the max heat-rate 
trajectory for the large ellipsled configuration using 
both LAURA and DPLR. Radiation transport calcula- 
tions were performed only with NEQAIR96 before the 
vehicle configuration evolved to the smaller, 2.88 m 
geometry. Centerline convective heating distributions 
(LAURA and DPLR results essentially identical) for 
each trajectory point are shown in Fig. 5 and global 
convective distributions at the peak heating point 
(t=180 s) along the trajectory are shown in Fig. 6. Peak 
convective and radiative heating rates (at the nose) for 
each trajectory point are shown in Fig. 7. The maxi- 
mum heating rates on the ellipsled were 3833 W/cm 
for convective heating and 1302 W/cm for radiative 
heating at t = 1 80 s. The highest percentage contribu- 
tion of radiation to the total heating environment was 
48% at t = 150 s. 



Figure 5: Centerline Convective Heating along Max 
Heat-Rate Trajectory for Large Ellipsled 



Figure 6: Global Convective Heating Distribution 
at Peak Heating for Large Ellipsled 



Figure 7: Nose Heating Rates along Max Heat-Rate 
Trajectory for Large Ellipsled 
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Uncertainties in Flow Field and Radiation Trans- 
port Methods 

Analysis of heating results and comparisons 
between the different flow field and radiation transport 
codes revealed several areas in which large uncertain- 
ties exist in the modeling of the high-energy aerocap- 
ture pass through Neptune’s atmosphere. Each of these 
areas will be discussed in the following subsections. In 
some cases differences were explored using a 0.3 m 
radius hemisphere geometry (the approximate nose 
radius of the large ellipsled) in place of the complex, 
three-dimensional ellipsled geometry, along with the 
simpler H 2 -He atmosphere. 

Flow Field Code-to Code Comparisons 

Code-to-code comparisons between LAURA and 
DPLR of convective heating rates resulted in agree- 
ment to within ±10% or less, as shown by sample com- 
parisons in Fig. 8 for the large ellipsled and Fig. 9 for 
the small ellipsled. As a result, where convective heat- 
ing rates for the ellipsled are shown throughout this 
paper, they are usually not identified as resulting from 
either code because the differences were generally very 
small. However, as will be shown subsequently, differ- 
ences in kinetic models initially implemented in the 
two codes which had little influence on convective 
heating rates had large effects on radiation computa- 
tions. 

Vibrational Nonequilibrium 

Stability problems were encountered using both 
LAURA and DPLR to compute the ellipsled flow fields 
in the H 2 -He-CH 4 Neptune atmosphere. These prob- 
lems were traced back to the modeling of vibrational 
non-equilibrium. For the trajectories under consider- 
ation, dissociation of the molecular species (H 2 , CH 4 , 
CH 3 , CH 2 , CH, C 2 ) behind the shock was very rapid. 
This dissociation left very small concentrations of mol- 
ecules for which vibrational equilibrium could be 
defined except in narrow regions at the shock and wall, 
which made the computational problem very stiff. 
Examples of this rapid dissociation are given by stagna- 
tion-line temperatures and H 2 mole fractions plotted in 
Fig. 10 for the 0.3 m hemisphere geometry at the peak 
heating points on the max heat load and max heat rate, 
|3 = 895 kg/m 2 trajectories. 



Figure 8: LAURA-DPLR Centerline Convective 
Heating Comparison for Large Ellipsled 



Figure 9: LAURA-DPLR Centerline Convective 
Heating Comparison for Small Ellipsled 



Figure 10: Hemisphere Stagnation-Line Tempera- 
ture and H 2 Mole Fractions 
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An additional uncertainty in vibrational non-equi- 
librium computations is the use of Park’s two-tempera- 
ture (TT V ) model 25 . This model was developed for use 
in air, where the primary diatomic species, N 2 and 0 2 , 
have similar vibrational characteristics. The applicabil- 
ity of this model to the Neptune environment is 
unproven. 

The flow field was found to be very nearly in 
vibrational equilibrium for points along the lower-alti- 
tude, max heat-rates trajectories, while vibrational 
non-equilibrium was present for points along the 
higher-altitude, max heat-load trajectories. In order to 
avoid the stability problems which occurred when 
vibrational non-equilibrium was allowed, all small 
ellipsled cases were computed with vibrational equilib- 
rium. This approximation had no appreciable affect on 
convective heating rates, but was recognized as a con- 
servative assumption in the computation of radiative 
heating rates. The conservatism results from the fact 
the equilibrium temperatures used to evaluate radiation 
transport were considerably higher near the shock than 
the vibrational temperatures which would have been 
used if vibrational non- equilibrium was allowed. 
These temperature differences are illustrated for the 0.3 
m hemisphere geometry (H 2 -He mixture) in Fig. 11. 

As an example of the differences resulting from 
the vibrational equilibrium assumption, the stagna- 
tion-point radiative heating levels computed using 
RADEQUIL with inputs from LAURA for the four 
cases shown in Fig. 11 were: 

Peak convective heat-load trajectory (t = 300 s): 

rs 

q ra d = 6084 W/cm , vibrational non-equilibrium 
Trad = 36,950 W/cm 2 , vibrational equilibrium 
Peak convective heat-rate trajectory (t = 184 s): 
Trad = 5192 W/cm 2 , vibrational non- equilibrium 
Trad = 7433 W/cm 2 , vibrational equilibrium 
As shown by the above values, the vibrational 
equilibrium assumption does reasonably well (for radi- 
ation calculations) at approximating the radiative heat- 
ing rates for cases near equilibrium, but for cases in 
non-equilibrium, this assumption produces extremely 
conservative results. 

The differences in radiative heating levels for the 
peak heat-load cases were clearly unacceptable; how- 
ever, the numerical stability problems discussed previ- 
ously prevented non-equilibrium solutions from being 
computed on the ellipsled geometry in a timely manner. 
Therefore, an approximate correction to the equilib- 
rium radiative heating results was made to obtain 
order-of-magnitude estimates for non-equilibrium lev- 
els. It was found that radiative heating levels approxi- 
mating those obtained from non-equilibrium flow field 
computations could be obtained from equilibrium flow 
field computations if the H 2 molecular transition contri- 


butions were neglected. The rationale for this approxi- 
mation was that H 2 produced the majority of the total 
radiation in the peak-load trajectory cases with equilib- 
rium modeling, but with non- equilibrium modeling, the 
H 2 would have been radiating at a much lower tempera- 
ture. Thus, subtracting the H 2 contribution from the 
equilibrium total radiative heating rate resulted in radi- 
ative heating estimates of the same order of magnitude 
as non-equilibrium predictions. For the hemisphere 
cases discussed above, the approximate (H 2 contribu- 
tion removed) equilibrium radiative heating rates were: 
Peak load trajectory (t = 300s): 

Trad = 5684 W/cm equilibrium with H 2 removed 
Peak rate trajectory (t = 184 s): 

Trad = 4797 W/cm 2 equilibrium with H 2 removed 



Figure 11: Hemisphere Stagnation Line Transla- 
tional and Vibrational Temperatures 

Radiation Code-to-Code Comparisons 

Comparison of LAURA-RADEQUIL and 
DPLR-NEQAIR96 flow field equilibrium radiative 
heating levels (with the H 2 contributions removed) 
were performed for several cases. 

For the hemisphere cases discussed in the previous 
section, the predicted levels were: 

Peak convective heat-load trajectory (t = 300 s): 
q rad = 5684 W/cm 2 LAURA-RADEQUIL 
Qrad = 1800 W/cm 2 DPLR-NEQAIR96 
Peak convective heat-rate trajectory (t = 184 s): 
qrad = 4797 W/cm 2 LAURA-RADEQUIL 
qrad = 4100 W/cm 2 DPLR-NEQAIR96 
For the small ellipsled, the peak (at the nose) radia- 
tive heating levels were: 

Peak convective heat-load trajectory (t = 300 s): 
q rad =8120 W/cm 2 LAURA-RADEQUIL 
q rad = 2200 W/cm 2 DPLR-NEQAIR96 
Peak convective heat-rate trajectory (t = 184 s): 
q rad = 5610 W/cm 2 LAURA-RADEQUIL 
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q iad = 4400 W/cm 2 DPLR-NEQAIR96 

From these numbers, the differences between the 
two sets of codes were found to be -25% for near-equi- 
librium conditions (along the peak heat-rate trajectory) 
but were up to -250% for non-equilibrium conditions 
(along the peak heat-load trajectory). 

As flow field code-to-code comparisons revealed 
only minor differences (when the same kinetic models 
were employed), the differences in radiative heating 
levels were attributed almost entirely to the radiation 
transport solvers. These differences remained unre- 
solved within the time-frame of this study, but several 
different assumptions in the radiation transport models 
of the two codes were noted: RADEQUIL includes 
more H atomic line transitions (Lyman-a,|3,y,6,e 
Balmer-a,|3,y,6, and Paschen-a,|3,y) than NEQAIR96 
(Lyman-a,p,y Balmer-a,p,y); NEQAIR96 includes 
line-by-line calculations of all radiation wavelengths, 
while RADEQUIL uses a smeared molecular band 
model; RADEQUIL includes more bound-free 
photo-ionization transitions (Lyman, Balmer, Paschen, 
Brackett and approximate integration thereafter to oo) 
than NEQAIR96 (Lyman and Balmer). 

Chemical Kinetics 

Comparisons of computations for the large 
ellipsled showed that the two flow field solvers pro- 
duced different results for shock stand-off distances 
and post-shock temperatures. While these differences 
had very little effect on the convective heating levels, 
they did lead to different predictions for radiative heat- 
ing levels. It was determined that the differences were 
due to the use of the Nelson kinetics in the LAURA 
computations and the Leibowitz kinetics in DPLR. 
These differences are illustrated in Figs. 12-13. The 
results shown are for the 0.3 m hemisphere using 
LAURA with both the Nelson or Leibowitz kinetics 
and an H 2 -He mixture. In Fig. 12, the stagnation line 
temperatures are shown along with the mole fraction of 
ionized hydrogen, the ionization rate of which was the 
main reason for the different temperatures predicted 
using the two kinetic models. As shown in Fig. 13, the 
different models led to only about a ±10% difference in 
convective heating. While the accuracy of both models 
for Neptune flow fields needs to be further explored, 
the Leibowitz kinetics led to much higher radiative 
heating predictions. Therefore, this more conservative 
model was employed in both codes for subsequent 
computations on the small ellipsled. 



Figure 12: Kinetic Model Effects on Hemisphere 
Stagnation-Line Temperature and H+ Mole Frac- 
tions 



Figure 13: Kinetic Model Effects on Hemisphere 
Convective and Radiative Heating Distributions 


Atmospheric Composition 

In this study, the baseline atmospheric composition 
of Neptune was assumed to 0.6246 for H 2 , 0.2909 for 
He and 0.0846 for CH 4 by mass. However, there is evi- 
dence to suggest a trace amount of N 2 in Neptune’s 
atmosphere. In order to determine the effects of com- 
position on the heating environment, 
DPLR-NEQAIR96 computations were performed for 
the large ellipsled at the peak heating point (t = 180 s) 
for 5 species (H 2 /He and products), 1 1 species 
(H 2 /He/CH 4 and products without C 2 or He±), and 19 
species (H 2 /He/CH 4 /N 2 and products with different N 2 
fractions) compositions. Centerline convective heating 
distributions and stagnation point convective and radia- 
tive heating rates for these cases are shown in Fig. 14. 
Convective heating levels were found to be relatively 
insensitive to composition, while radiative heating lev- 


7 


els were sensitive to the presence of N 2 through the for- 
mation of radiating CN. 



Figure 14: Effects of Atmospheric Composition on 
Convective and Radiative Heating 


Radiation-Flow Field Coupling 

In this study, radiation and flow field computations 
were uncoupled. However, previous studies (e.g. 26) 
have shown that coupling of computations (feeding the 
radiation transport results back into the flow field code) 
can result in significant reductions to the predicted radi- 
ative heating levels through non-adiabatic radiative 
cooling of the flow field. The current uncoupled 
approach is thus recognized to yield conservative 
results. 

2.88 m, P = 895 kg/m 2 Vehicle 

Convective and radiative heat transfer computa- 
tions were performed for seven points (Table 2) along 
the max heat-load trajectory and at the peak heating 
point on the max heat-rate trajectory. As per the dis- 
cussion in the previous section, the flow field was mod- 
eled as being in vibrational equilibrium, while 
approximate radiative heating rates were computed by 
neglecting the H 2 contribution. Centerline convective 
heating distributions for each trajectory point are 
shown in Figs. 15 and 16. Peak (at the nose) convec- 
tive and radiative heating rates for each trajectory point 
are shown in Fig. 17. The maximum convective heat- 
ing rates were 2575 W/cm on the max-load trajectory 
and 7915 W/cm on the max-rate trajectory. Consider- 
able differences were again observed in the radiative 
heat transfer rates from RADEQUIL and NEQAIR96, 
but both sets of results showed that radiative heating 
was of the same order-of-magnitude as convective 
heating. At the peak heating point on the 
max-heat-load trajectory, DPLR-NEQAIR96 predicted 
radiation heating levels at -45% of the total heating 


while LAURA-RADEQUIL predicted radiative levels 
at -80% of the total heating level. 



Figure 15: Centerline Convective Heating along 
Max Heat-Load Trajectory for Small Ellipsled 



Figure 16: Centerline Convective Heating along 
Max Heat-Rate Trajectory for Small Ellipsled 
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TPS DEVELOPMENT 

Approximate heating rate and integrated heat-load 
environments were generated from these results for use 
in TPS material selection and thickness sizing. The 
convective and radiative rates from the computations 
were used as anchor points from which to scale the 
heating time-history outputs from POST for a 1-m 
hemisphere to the ellipsled geometry. Three environ- 
ments were generated: a “Low” environment based on 
the convective heating plus one-half of the radiative 
heating rate to account for radiative cooling effects, a 
“Reference” environment based on the total (convec- 
tive plus radiative) heating rates generated herein, and a 
“High” environment based on twice the total heating 
rates to account for turbulent heating augmentation. 
TPS development based on these environments is dis- 
cussed in the companion paper by Laub 5 . 

SUMMARY 

Preliminary convective and radiative heating envi- 
ronments for a Neptune aerocapture mission have been 
computed. Environments were generated both for a 
large 5.50 m ellipsled and a small 2.88 m ellipsled. 
Radiative heating constituted up to 80% of the total 
heating along the trajectories studied. 

Because of the expected computational difficul- 
ties for this high-velocity aerocapture mission in Nep- 
tune’s H 2 -He-CH 4 atmosphere, heating environments 
were generated in tandem using LAURA with RADE- 
QUIL and DPLR with NEQAIR96 to compute the flow 
field and radiation transport properties. This approach 
was designed to reduce uncertainties and to identify 
areas in which further research and development of 
numerical models and tools will be required in order to 
provide higher confidence in analyses for this class of 
mission. 

The computations were found to agree well for 
flow field properties and convective heating distribu- 
tions (when the same kinetic models were employed), 
but several sources of large uncertainty were identified 
in the computation of radiative heating. 

Kinetic modeling of reactions in the H 2 -He-CH 4 
Neptune atmosphere was one of the problem areas 
identified. The use of different reaction sets cited in the 
literature produced large differences in post-shock 
stand-off distance and temperatures, which led to large 
differences in predicted radiative heating rates. 

Vibrational non-equilibrium modeling also pre- 
sented difficulties. Numerical stability could not be 
achieved for several cases when a two-temperature, 
vibrational non-equilibrium model was employed. In 


order to resolve this issue, vibrational equilibrium was 
imposed, which was shown to lead to over-prediction 
of radiative heating levels. Furthermore, even if stabil- 
ity had been achieved, the use of a two-temperature 
model developed for Earth’s N 2 -0 2 atmosphere is 
unproven in Neptune’s H 2 -He-CH 4 atmosphere. 

It was also found that large differences existed in 
radiative heating rates produced by the two radiation 
transport codes. While several differences in the physi- 
cal models incorporated in these codes were identified, 
the specific reasons for these differences were not iden- 
tified. 

Within the limits of this short-term study, prelimi- 
nary convective-radiative heating environments for 
thermal protection system sizing were generated which 
were sufficient to support a moderate-fidelity vehicle 
design for a Neptune aerocapture mission. However, in 
order to complete a high-fidelity design, further devel- 
opment of computational tools and methods for the 
Neptune environment will be required. 

ACKNOWLEDGEMENTS 

This research was supported by the In- Space Pro- 
pulsion Program Office at the NASA Marshall Space 
Flight Center. Dinesh Prabhu was supported by NASA 
Ames Research Center grant NAS2-99092 to ELORET 
Corp. Naruhisa Takashima was supported by NASA 
Langley Research Center through contract 
NAS 1-00 135 to AMA, Inc. Kenneth Sutton was sup- 
ported through NASA Langley contract NCC 1-02043 
to the National Institute of Aerospace. 

REFERENCES 

Lockwood, M. K., “Neptune Aerocapture Systems 
Analysis”, AIAA Paper 2004-4951, Aug. 2004. 
2 Edquist, K. T, Prabhu, R., Hoffman, D., and Rea, J., 
“Configuration, Aerodynamics, and Stability Analysis 
for a Neptune Aerocapture Orbiter,” AIAA Paper 
2004-4953, Aug. 2004. 

3 Starr, B., Powell, R., Westhelle, C., and Masciarelli, 
J., “Aerocapture Performance Analysis for a Nep- 
tune-Triton Exploration Mission,” AIAA Paper 
2004-4955, Aug. 2004. 

4 Justus, C, Duvall, A., and Keller, V., “Atmospheric 
Models for Aerocapture Systems Study,” AIAA Paper 
2004-4952, Aug. 2004. 

5 Laub, B. and Chen, Y., “Challenges for Neptune 
Aerocapture,” AIAA Paper 2004-5178, Aug. 2004. 
6 Powell, R.W., Striepe, S.A., Desai, P.N., Queen, E.M.; 
Tartabini, P.V., Brauer, G.L., Cornick, D.E., Olson, 
D.W., Petersen, F.M., Stevenson, R., Engel, M.C.; 
Marsh, S.M., “Program to Optimize Simulated Trajec- 
tories (POST II), Vol. II Utilization Manual,” Version 


9 



1.1. l.G, May 2000, NASA Langley Research Center, 
Hampton, VA. 

7 Sutton, K. and Graves, R. A., “A General Stagna- 
tion-Point Convective-Heating Equation for Arbitrary 
Gas Mixtures,” NASA TR-R-376, Nov. 1971. 

8 Gnoffo, P. A., “An Upwind-Biased, Point-Implicit 
Algorithm for Viscous, Compressible Perfect-Gas 
Flows,” NASA TP-2953, Feb. 1990. 

9 Cheatwood, F. M., and Gnoffo, P. A., “User’s Manual 
for the Fangley Aerothermodynamic Upwind Relax- 
ation Algorithm (FAURA),” NASA TM 4674, April, 
1996. 

10 Wright, M.J., Candler, G.V., and Bose, D., 
“Data-Parallel Line Relaxation Method for the 
Navier-Stokes Equations,” AIAA Journal , Vol. 36, No. 
9, 1998, pp. 1603-1609. 

n Roe, P. L., “Approximate Riemann Solvers, Parame- 
ter Vectors and Difference Schemes,” Journal of Com- 
putational Physics , Vol. 43, No. 2, 1981, pp. 357-372. 

1 9 

Harten, A., “High Resolution Schemes for Hyper- 
bolic Conservation Laws,” Journal of Computational 
Physics , Vol. 49, No. 3, 1983, pp. 357-393. 

13 Yee, H. C., “On Symmetric and Upwind TVD 
Schemes,” NASA TM 88325, 1990. 

14 MacCormack, R.W. and Candler, G.V., “The Solu- 
tion of the Navier-Stokes Equations Using 
Gauss- Seidel Line Relaxation,” Computers and Fluids , 
Vol. 17, No. 1, 1989, pp. 135-150. 

15 Yee, H.C., “A Class of High-Resolution Explicit and 
Implicit Shock Capturing Methods,” NASA TM 
101088, Feb. 1989. 

16 Takashima, N., Hollis, B., Olejniczak, J., Wright, M., 
and Sutton, K., “Preliminary Aerothermodynamic s of 
Titan Aerocapture Aeroshell,” AIAA Paper No. 
2003-4952, July 2003. 

17 01ejniczak, J., Wright, M., Prabhu, D., Takashima, 
N., Hollis, B., Zoby, E. V., and Sutton, K. “An Analy- 
sis of the Radiative Heating Environment for Aerocap- 
ture at Titan,” AIAA Paper 2003-4953, July 2003. 
18 Nelson, H. F., Park, C., and Whiting, E. E., “Titan 
Atmospheric Composition by Hypervelocity 
Shock-Layer Analysis,” Journal of Thermophysics and 
Heat Transfer , Vol. 5, No. 2, April-June 1991, pp. 
157-165. 

19 Park, C., “Radiation Enhancement by Nonequilib- 
rium During Flight Through the Titan Atmosphere,” 
AIAA Paper 1982-0878, June 1982. 

20 Leibowitz, L. P., and Kuo, T., “Ionizational Nonequi- 
librium Heating During Outer Planetary Entries,” AIAA 
Journal , Vol. 14, No. 9, Sept. 1976, pp. 1324-1329. 
21 McBride, B. J., Zehe, M. J., and Gordon, S., “NASA 
Glenn Coefficients for Calculating Thermodynamic 
Properties of Individual Species,” NASA TP 
2002-211556, Sept. 2002. 


22 Nicolet, W. E., “User’s Manual for the Generalized 
Radiation Transfer Code (RAD/EQUIL),” 
NASA-CR-1 16353, Oct. 1969. 

23 Nicolet, W. E., “Advanced Methods for Calculating 
Radiation Transport in Ablation-Product Contaminated 
Boundary Layers,” NASA-CR-1 656, Sept. 1970. 
24 Whiting, E. E., Yen, L., Arnold, J. O. and Paterson, J. 
A., “NEQAIR96, Nonequilibrium and Equilibrium 
Radiative Transport and Spectra Program: User’s Man- 
ual,” NASA RP-1389, Dec. 1996. 

25 Park, C., “Assessment of Two Temperature Kinetic 
Model for Ionizing Air,” AIAA Paper 1987-1574, June 
1987. 

26 Wright, M. J., Bose, D., and Olejniczak, J., “The 
Impact of Flowfield-Radiation Coupling on Aeroheat- 
ing for Titan Aerocapture,” AIAA Paper 2004-0484, 
Jan. 2004. 


10 



Table 1: Free Stream Conditions for Large Ellipsled Trajectory Points 


Trajectory 

Time 

(s) 

Altitude 

(m) 

Density 

(kg/m3) 

Temperature 

(K) 

Velocity 

(m/s) 

Max convective heat rate 

150 

207,090 

1.319E-05 

132.41 

31,450 

Max convective heat rate 

170 

148,079 

8.392E-05 

106.95 

30,534 

Max convective heat rate 

180 

132,186 

1.450E-04 

103.16 

29,243 

Max convective heat rate 

190 

130,444 

1.538E-04 

102.95 

27,670 

Max convective heat rate 

210 

160,550 

5.397E-05 

114.43 

25,777 


Table 2: 

Free Stream Conditions for Small Ellipsled Trajectory Points 


Trajectory 

Time 

(s) 

Altitude 

(m) 

Density 

(kg/m3) 

Temperature 

(K) 

Velocity 

(m/s) 

Max convective heat rate 

184 

107,702 

3.513E-04 

96.95 

29,158 

Max convective heat load 

160 

389,739 

3.670E-06 

187.50 

31,524 

Max convective heat load 

170 

370,114 

5.149E-06 

186.60 

31,506 

Max convective heat load 

180 

352,568 

7.016E-06 

185.70 

31,478 

Max convective heat load 

200 

323,363 

1.239E-05 

183.70 

31,374 

Max convective heat load 

300 

261,092 

3.709E-05 

180.01 

30,049 

Max convective heat load 

450 

260,262 

3.768E-05 

179.93 

27.309 

Max convective heat load 

600 

299,445 

1.822E-05 

182.67 

25,581 
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Table 3: Nelson-Park Kinetic Model 


# 

n -(9/TJ 
k f = CT" e 

C 

(cc/mol/s) 

n 

e 

(K) 

T a 

(K) 

Ref. 

1 

ch 4 + m^ch 3 + h + m 

2.25x1 0 27 

-1.87 

52,900 

(TT V ) 0 ' 3 

18-19 

2 

ch 3 + m^ch 2 + h + m 

2.25xl0 27 

-1.87 

54,470 

(TT V ) 05 

18-19 

3 

ch 2 + m^ch + h + m 

2.25xl0 27 

-1.87 

50,590 

(TT V ) 0 ' 5 

18-19 

4 

CH + M^C + H + M 

1.13xl0 19 

-1.00 

40,193 

(TT V ) 0 ' 5 

18-19 

5 

c 2 + m^c + c + m 

9.68xl0 22 

-2.00 

71,000 

(TT V ) 0 ' 5 

18-19 

6 

h 2 + m^h + h + m 

1.47xl0 19 

-1.23 

51,950 

(TT V ) 05 

18-19 

8 

H, + C -* CH + H 

1.80xl0 14 

0.00 

11,490 

T 

18-19 

9 

L + . + 

C + e — ^ C + e + e 

3.90xl0 33 

-3.78 

130,000 

T 

A V 

18-19 

10 

H + e H + + e + e + 

5.90xl0 37 

-4.00 

157,800 

T 

A v 

18-19 

11 

He + e -> He + + e + e + 

1.33xl0 13 

0.50 

286,160 

T 

A V 

20 


Table 4: Leibowitz-Nelson-Park Kinetic Model 


# 

n -(9/T ) 
k f = CT^ e 

c 

(cc/mol/s) 

n 

e 

(K) 

T a 

(K) 

Ref. 

1 

ch 4 + m^ch 3 + h + m 

2.25x1 0 27 

-1.87 

52,900 

(TT v ) o;> 

18-19 

2 

CH 3 + M -»• CH 2 + H + M 

2.25xl0 27 

-1.87 

54,470 

(TT V ) 0 ' 5 

18-19 

3 

ch 2 + m^ch + h + m 

2.25xl0 27 

-1.87 

50,590 

(TT V ) 05 

18-19 

4 

CH + M-»C + H + M 

1.13xl0 19 

-1.00 

40,193 

(TT V ) 0 ' 5 

18-19 

5 

c 2 + m^c + c + m 

9.68xl0 22 

-2.00 

71,000 

(TT V ) 05 

18-19 

6 

h 2 + m^h + h + m 

1.04xl0 19 

-1.00 

51,950 

(TT V ) 0 ' 5 

20 

8 

h 2 + c^ch + h 

1.80xl0 14 

0.00 

11,490 

T 

18-19 

9 

C + e — ^ C + e + e 

3.90xl0 33 

-3.78 

130,000 

T 

A V 

18-19 

10 

H + e— >H -he -he 

2.28xl0 13 

0.50 

157,800 

T 

A V 

20 

11 

H + e' -*« H + + e‘ + e + 

4.11xl0 13 

0.50 

116,100 

T 

A V 

20 

12 

H + H H + + e + H 

6.17xl0 10 

0.50 

116,100 

T 

20 

13 

H + He H + + e + He 

4.88xl0 10 

0.50 

116,100 

T 

20 

14 

He + e -> He + + e + e + 

1.33xl0 13 

0.50 

286,160 

T 

A v 

20 


Note: Equations 10 and 1 1 in Table 4 are one-step and two-step electron-impact ionization reactions for H 
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